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Abstract 

We investigate a model of brane cosmology to find a unified description of the radiation- 
matter-dark energy universe. It is of the interacting holographic dark energy with a bulk- 
holographic matter x- This is a five-dimensional cold dark matter, which plays a role of 
radiation on the brane. Using the effective equations of state oo'^^ instead of the native 
equations of state cua, we show that this model cannot accommodate any transition from 
the dark energy with u'^ > —1 to the phantom regime u'^ < —1. Furthermore, the case 
of interaction between cold dark matter and five dimensional cold dark matter is consid- 
ered for completeness. Here we find that the redshift of matter-radiation equality z^^ is 
the same order as = 2.4 x lO'^Q^^h^. Finally, we obtain a general decay rate F which is 
suitable for describing all interactions including the interaction between holographic dark 
energy and cold dark matter. 



*e-mail address: ysniyung@inje.ac.kr 



1 Introduction 



Recent observations from Supernova (SN la) P and large scale structure [2] imply that our 
universe is accelerating. Also cosmic microwave background observations [31 H] provide an 
evidence for the present acceleration. A combined analysis of cosmological observations 
shows that the present universe consists of 70% dark energy and 30% dust matter including 
cold dark matter (CDM) and baryons. 

Although there exist a number of dark energy models, a promising candidate is the 
cosmological constant. However, one has the two famous cosmological constant problems: 
the fine-tuning and coincidence problems. In order to solve the first problem, we may 
introduce a dynamical cosmological constant model inspired by the holographic principle. 
The authors in [5] showed that in quantum field theory, the UV cutoff A could be related 
to the IR cutoff La due to the limit set by introducing a black hole (the effects of gravity). 
In other words, if pa = is the vacuum energy density caused by the UV cutoff, the total 
energy of system with the size La should not exceed the mass of the black hole with the 
same size La: LaPa < SM^La. If the largest cutoff La is chosen to be the one saturating 
this inequality, the holographic energy density (HDE) is given by PA = 3c^M^/87rLl with 
a constant c. The lower limit of c is protected as c > 1 by the entropy bound. Here 
we regard pa as a dynamical cosmological constant. Taking the IR cutoff as the size 
of the present universe (La = ^/H), the resulting energy is close to the present dark 
energy [6]. However, this approach with L/^ = 1/ H is not fully satisfied because it fails to 
recover the equation of state (EoS) for the dark energy-dominated universe [7]. Further 
studies in [SI El [IHl E] have shown that choosing the future event horizon as the IR cutoff 
determines an accelerating universe with the native EoS cua = — 1/3(1 + cilnpA/cilna) = 
-l/3-2x/ni/3c. 

Also if the interaction is turned on, the coincidence problem could be resolved [12] . 
The interacting dark energy models provided a new direction to understand the dark 
energy [T311I11II3]. The authors in [TB] introduced an interacting holographic dark energy 
model where an interaction exists between HDE and CDM. They derived the phantom- 
phase of C(Ja < — 1 using uj\. However, it turned out that the interacting holographic dark 
energy model could not describe a phantom regime when using the effective equation of 
state ci^A^ [XT]. More recently, it was shown that for non-fiat universe of 7^ [TSl [TU] , 
the interacting holographic dark energy model could not describe a phantom regime of 
u]^ < — 1 [Sn]. In Ref . [2T] . the authors discussed the cosmological dynamics of interacting 
holographic dark energy model using the phase-space variables. A key of this system is an 
interaction between two matters. Their contents are changing due to energy transfer from 
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HDE to CDM until the two components are comparable. If there exists a source/sink in 
the right-hand side of the continuity equation, we must be careful to define its EoS. In 
this case, the effective EoS is the only candidate to represent the state of the mixture of 
two components arisen from decaying of HDE into CDM. This is clearly different from 
the non-interacting case which can be described by the native EoS completely. 

On the other hand, if the brane cosmology is introduced, one could have interesting 
interaction between bulk and brane matters. In the low energy limit, the brane cosmology 
reduces to the Friedmann-Robertson- Walker (FRW) form with a bulk-holographic matter 
X- This is just a five dimensional cold dark matter(5DCDM) which play a role of a four- 
dimensional radiation when using the effective EoS approach. Then a unified description 
of radiation-matter-dark energy universe could be performed within the brane cosmology. 
Here we obtain two kinds of interaction: HDE-5DCDM and CDM-5DCDM. The first 
interaction may be allowed because one may allow the interaction of HDE with radiation. 
However, the latter seems not to be permitted because we assume that the CDM is not a 
source of radiation and it does not interact with the radiation. However, we suggest that 
the two interactions are possible to occur within the brane cosmology. 

Concerning the brane-bulk interaction, there were contradictions: if one uses the effec- 
tive EoS of w^f 0, a transition occurs between w^f > —1 and w'^ < —1 [22]. On the other 
hand, using w^^, it was shown that such a transition does not occur In this work, we 
wish to address this issue again. We solve three coupled differential equations for density 
parameters numerically by assuming three interactions between them. Furthermore, 
we introduce three types of the decay rate F to find the dark energy-dominated evolution 
on the brane. We confirm that any phantom-phase is not found on the brane. 

2 Brane-bulk interaction model 

Generalization of the Randall- Sundrum scenario [24J in cosmology considers the AdSs 
geometry containing the bulk cosmological constant A, but explores arbitrary energy 
densities on the brane and in the bulk. The Binetruy-Deffayet-Langlois (BDL) approach 
is a genuine extension of the Kaluza-Klein cosmology to account for the local distribution 
on the brane [23]. In this case, the location of the brane is fixed with respect to the bulk 
direction. This approach is useful for describing the cosmological evolution of the brane 

^The authors in [5^ use a different definition w'^ = — f — g 'rfin^ from our definition w'^^ . Here 
5H^ — / Hq — r^m/a^ accounts for all terms in the Friedmann equation not related to the brane matter 
fim- The = — f crossing is achieved by considering the brane-bulk interaction without specifying 
dark energy as holographic dark energy. 
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when a brane-bulk interaction exists. Hence, we follow the BDL brane cosmology. We 
introduce the gaussian-normal bulk metric for (1 + 3 + l)-dimensional spacetime 



d^BDL = -c {t,y)dt +a {t,y)%jdx'dx^ + b {t,y)dy , (1) 

where jij is the metric of a three-dimensional space with a constant curvature of 6k. Let 
us express the bulk Einstein equation Gmn = -2M^Tmn in terms of the BDL metric §. 
We introduce a (1 + 3)-dimensional brane located at ?/ = 0. For simplicity, we choose 
the total stress-energy tensor = diag(— A, —A, —A, —A, —A) + n + u- Here 

A is the bulk cosmological constant and the bulk stress-energy tensor T*^ n from 
is not needed to have a specific form initially. If T* = 0, it is obvious that there is 
no brane-bulk interaction. The brane stress-energy tensor from including the brane 
tension a and the energy density p is assumed to take the form 

u = -^dmg{-p - a,p - a,p - a,p - a,0). (2) 

We are interested in solving the Einstein equations at the location of the brane. Initially 
we indicate by the subscript "0" for the value of various quantities on the brane. Also it 
is convenient to choose the gaussian-normal gauge with Bq = 1 and the temporal gauge 
with Co = 1 on the brane. We obtain from G^y = jM^'^oy, 

p + 3%(l + a;) = -2T%. (3) 
ao 

Here we assumed an equation of state p = up on the brane. 

On the other hand, the average part of yy-component equation is given by 

Then, we rewrite Eq.dl]) in the following equivalent form by introducing the two bulk- 
holographic energy densities x ^i-nd (p: 

^ + ^HoX = ^Ap + ^)f\, (6) 



36M6 
1 

3M3' 



+ 4Hoct> = -^Hofy y, (7) 



^Our action is given by S'g = / d^Xy/^{APR - A + C'g^^j + J d'^xyf^O^"^ with AP = l/ieTrGs 
1/2^1 and rjj'"* = -(cr + p) [ig. 
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with Hq = do/oo. In the case of p = p = and = x = 0, one finds the Randall- Sundrum 
vacuum state [21]. We choose the cosmological constant A = — cr^/12M^ = — 12M^/£^ 
with the brane tension a = 12M^/i to have a critical brane. Hence the cosmological evo- 
lution will be determined by four initial parameters (pi, aoi, Xi, (pi) instead of two (pi, aoi) 
in the FRW universe. This is so because the generalized Friedmann equation ([5]) is not 
a first integral of the Einstein equation. It is mainly due to the energy exchange T* y 
between the brane and bulk. In the case of = and y = Ap > with T"^ y = 0, 
one finds a mirage-radiation term x ~ (1 — e~^*''^)/ao for an energy outfiow from the 
brane [27] . It is a cosmological model that the real matter on the brane decays into the 
extra dimension. Also for = 0,T* ^ ~ — ^ and ^ = 0, it is shown that the energy 
influx from the bulk generates a cosmological acceleration on the brane with the acceler- 
ation parameter Q = 7^f^ = l — f? where < g < 4 [28]. However, in general, it will be 
a formidable task to solve Eqs.(l5])-(I7j) with Eq.Q because it gives rise to a complicated 
dynamics between the brane and the bulk. In Ref.[22], they used T* ^ oc if a" to derive 
the super-acceleration using u^. 

For our purpose, let us imagine a brane universe made of CDM pm with cUm = 0, but 
obeying the holographic principle. In addition, we propose that the holographic energy 
density p\ exists with its native EoS cua > — 1 on the brane. If one assumes a form of the 
interaction T with (j) = y = 0, their continuity equations take the simple formj^l 

P + 3H{1 + Uj)p = -T, p = PA + Pra (8) 

X + ^Hx = T (9) 

and the generalized Friedmann equation ([5]) on the critical brane leads to 

87r r 1 A; 



3M2 



p + x --• (10) 



Now we consider the case of decaying from HDE to 5DCDM with T = FpA, while the 
CDM is conserved by choosing 

PA + 3H{1 + u;a)pa = -T, (11) 
p^ + 3Hp^ = 0. (12) 

This decaying process impacts their equations of state and particularly, it induces the 
effective EoS for the 5DCDM. Interestingly, an accelerating phase could arise from a large 



■^Hereafter, we focus on the brane. Hence we use the notation without the subscript "0" and T = 2T" 



and X — (72M^/(t)x, and a/77M^ = 1/6M^£ = 87r/3Mp. Also we concentrate on the fow-energy region 
of p <C cr and thus p^-term in Eq.© is neghgible. 
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effective non-equilibrium pressure 11^ defined as = —Tp/^/3H{= Ha). Then the two 
equations (fTT|) and © are translated into those of the two dissipatively imperfect fluids 



PA + 3H l + tUA + — 

L oil 



Pa = Pa + 3iJ (1 + u;a)pa + IIa 



0, 



X + 3H 



1 + 

3 X 



x-x + 3i^i(i + 3)x-n. 



0. 



(13) 
(14) 



The positivity of IIa > shows a decaying of HDE via the cosmic frictional force, while 
Uy. < induces a production of the mixture via the cosmic anti-frictional force simultane- 
ously [29l|30]. This is a sort of the vacuum decay process to generate a particle production 
within the two-fluid model j3_l j . As a result, a mixture of two components will be created. 
When turning on the interaction term, from Eqs. fll3p and fll4p . we read off their effective 
equations of state as 

cff , r g 1 Pa r 
uT =uja + — , 3 



(15) 



3H' 3 X 

Hence it is clear that the 5DCDM x plays a role of radiation on the brane, if there is no 
interaction. Introducing the density parameters defined by fli = p-J Pc as 

SyrpiT, „ SttjOa „ k 



3M2/72' 



»7rpA 
3M2if2^ 



3M2if2' 



(16) 



we can rewrite the Friedmann equation (iTOl) as a simplified form 



+ l^A ~l~ ^x — ^ ~^ ^k- 

Hereafter we use this relation instead of Eq. lfTOl) . 

For the non-flat universe of A; 7^ 0, we introduce the future event horizon L\ 
a^Fnit) = a^Fuit) with 

eFH(t) = 

Here the comoving horizon size is given by 

dr 

'0 Vl — kr"^ 



(17) 



R 



FH 



dt 
t a 



smn 



\k\r{t) 



(19) 



where leads to ^fh^(^) = sin V(t), ^fh°(^) = "^i^)^ and ^{t) = sinh V(t). For our 
purpose, we use a comoving radial coordinate r(t). 



r(t) 



smn 



(20) 
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La = ar(t) is a useful length scale for the non-flat universe [IB]. Its derivative with respect 
to time t leads to 

c 

La = HLa + af = —7= — cosny, (21) 

where cosn?/ = cosy, y, coshy for = 1, 0, —1 with y = \^Rfu/cl- Hereafter we consider 
three classes of interactions: HDE-CDM, HDE-5DCDM, and CDM-5DCDM. Using the 
definition of p\ = ^^^^2'' and (fT3|l . one finds the equation of state for HDE 



3 3c 

Here we can read off the effective EoS for HDE as 



1 2 /O 

Pa + SH\i ^-^cosnyl Pa = 0. (22) 



1 2JnA{x) 

u;f{x) = --- ^3^^ cosny (23) 

with X = Ina. At the first sight, the above effective EoS seems not to be relevant to the 
interaction, but it depends on the decay rate T through f^A- 



3 Unified picture for interactions 

For the interaction between HDE and CDM on the brane, we assume to have 

Pra + 3Hp^ = f , (24) 

PA + 3H{1 + uoa)pk = -f, (25) 
X + AHx = 0, (26) 

where T is chosen as T = FpA for decaying from HDE to CDM, while T = — Fpm for 
decaying from CDM to HDE. This case is not realized by the brane cosmology because the 
interaction T is effective on the brane. Hence there is no brane-bulk interaction (T = 0). 
However, we include this type of interaction for completeness. 

In the case of interaction between HDE and 5DCDM, their continuity equations are 
given by 

p^ + 3Hp^ = 0, (27) 
PA + 3if (1 + u;a)pa = -T, (28) 
X + AHx = T, (29) 

Here T is chosen as T = FpA for decaying from HDE to 5DCDM, whereas T = —Fx for 
decaying from 5DCDM to HDE. 
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Figure 1: (color online) Graph for the noninteracting case. For = and c = 1, k = 1 
evolution of fl\ (green), (red), and (blue) and the equations of state, uj\ (cyan) 
and u'^ = 1/3 (yellow) with = 0. Here a; = Ina moves backward direction (— ) or 
forward direction (+), starting at the present time a: = 0(ao = 1). 



Finally, the case of interaction between CDM and 5DCDM takes the form 

PA + 31^(1 + uja)pa = 0, 



(30) 
(31) 
(32) 

Tx for 



where T is chosen as T = Fpm for decaying from CDM to 5DCDM, while T : 
decaying from 5DCDM to CDM. 

By choosing appropriate effective equations of state, the above equations for all three 
cases can be unified as follows: 



p^ + 3H{1 + <ff)p^ = 0, 
PA + 3if(l + Ujf)pA = 0, 
X + SH{l + ujf)x = 0, 



(33) 
(34) 
(35) 



All effective EoS are summarized on the Table [TJ However, u'j^ is the same for all cases 
as is given by Eq. (l23|) . Here we choose three types for the decay rate F with 



(1) -type :F = 3^2(1 + ^), 



0.2: 

(36) 
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Table 1: Summary of effective equations of state and related information. Here I, A, m, 
and X represent interaction, HDE, CDM, and 5DCDM, respectively. The redshift factor 
Zeq is determined from the relation x — — ln(l + 2;) when Om 
available" . Finally, yes (no) represent the status of evolution. 
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Fig. 4b 
Fig. 4d 
Fig. 4f 
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(3) - type : T = 3Hb\l + 



(37) 
(38) 



(l)-type is known as a conventional form for the interaction between HDE and CDM. 
However, choosing this form leads to an unwanted evolution and thus we have to introduce 
another interaction (2)-type for the evolution of the dark energy-dominated universe. 
Finally, (3)-type is chosen because (2)-type is not suitable for describing the interaction 
between CDM and 5DCDM. Another types are found in Ref.|12]. 

In order to obtain differential equations for density parameters, Q^, and which 
govern evolution of the universe, we introduce 



fh 

Pc 



fii,i = m, A,x. 



(39) 



Differentiating Ri with respect to cosmic time t and then using appropriate definitions, 
we obtain three equations 



2 + (1 + 3<")fi^ + (1 + 3cuf )fiA + (1 + 3cuf )fi 
2 + (1 + ?,uof,)n„, + (1 + 3cuf )fiA + (1 + 3c^f )a 



3fi„,(l 



.cS\ 



(40) 



3r]A(i + 0, (41) 



2 + (1 + 3<^)1]„, + (1 + 3o;f )fiA + (1 + 3cu^")r]^ - ?>^^{l + cu^"), (42) 



.cff 



,cff\ 



,cff\ 



where ' is the differentiation with respect to a; = In a. These equations come from the 
first and second Friedmann equations combined with their continuity equations. In order 
to obtain solution, we have to solve the above coupled equations numerically by consid- 
ering the initial condition at present time^: ^'{^\x=q > 0, ^a — ^■'^'^i^k=i — 0.01, = 
0.28,nl = 0.01. 

The noninteracting case with 6^ = is depicted at Fig. 1, which shows the standard 
evolution for the HDE. Here the effective EoS reduces to the native EoS because of 
the absence of interactions except u'^ = 1/3 for 5DCDM x- Each matter satisfies its 
continuity equation. We find a sequence of dominance in the evolution of the universe: 
radiation^ CDM — i>dark energy. The redshift factor Zeq = 27.1 is determined from the 
relation of x = — ln(l + z) when f2m = Q^. 

Figs. 2a-f show the evolution for the interaction between HDE and CDM on the 
brane. The left column of HDE^CDM was already known but the right column shows 
new results. These all indicate evolutions for dark energy-dominated universe except 



''Here we use the data from the combination of WMAP3 plus the HST key project constraint on 

Ho m- 
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the case of CDM^HDE with the decay rate T = 3Hb'^{l + n^/Vl^) ^ Ej. This case 
provides a negative density parameter fl-^ < for the future evolution and thus induces 
the unwanted case of > 1. 

Figs. 3a-f indicate the evolution for the interaction between HDE and 5DCDM. This 
corresponds to the case of interaction between HDE and radiation on the brane. The 
left column is for HDE— > 5DCDM. An evolution for dark energy-dominated universe is 
possible for only the decay rate of (3)-type: F = 3Hb'^{l + Q\/Q^)Q\Q^. The right 
column is for 5DCDM— i>HDE. Here evolutions come out when choosing (2)and (3)-type. 
All forward evolutions are possible, whereas backward evolutions are not possible for 

(1) -type and HDE^ 5DCDM with (2)-type. 

Figs. 4a-f show the evolution for the interaction between CDM and 5DCDM. This 
corresponds to the case of interaction between CDM and radiation on the brane. The left 
column is for CDM-^ 5DCDM and the right column is for SDCDM^CDM. An evolution 
for dark energy-dominated universe is possible for only the decay rate of (3)-type. All 
backward evolutions seem not to be possible for (1) and (2)-types. Especially, we find the 
unwanted backward evolution of Q^a < 0, > 1 for the 5DCDM-^CDM with (1) and 

(2) -types. In this sense, (3)-type is considered as the general form of decay rate F. 

4 Discussions 

We investigate a unified description of radiation-matter-dark energy universe within the 
brane cosmology. It is confirmed that there is no phantom phase from brane-bulk in- 
teractions (HDE-5DCDM, CDM-5DCDM) and interaction on the brane (HDE-CDM) 
when using uj^. Thus our results favors Setare's case [23] but disfavors Cai-Gong- Wang's 
case [22]. This arises mainly because we used a different definition for the effective EoS 
uj]^ from Cai-Gong- Wang's case of as well as the HDE as dark energy. Recently, the 
authors in [32] showed that the interacting holographic dark energy with CDM may lead 
to the phantom phase using the native EoS ua. Also the authors in |33j showed that the 
brane-bulk interaction without the holographic dark energy accommodates the u = —1 
crossing when using u^. Hence, the issue is to choose an appropriate EoS for describing 
the dark energy universe. 

Also, we obtain an additional information from the unified picture of interactions. We 
suggest a sequence of the evolution: radiation-dominated universe — > matter-dominated 
universe —>■ dark energy-dominated universe. The 5DCDM plays the same role as a 
radiation on the brane. As is shown Fig.l and Table 1 , we have z^q = 27.1 which is 
not close to = 2.4 x Q^h"^ ~ 4.8 x 10'^ if there is no interaction. Interestingly, as 
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Figure 2: (color online) Six graphs for the interaction between HDE and CDM. For 
6^ = 0.2 and c = 1, k = 1 evolution of (green), (red), and (blue) and the 
effective equations of state, u"^ (cyan) and cu'^ = 1/3 (yellow) with a;^(pink). The left 
column is for HDE— >CDM and the right one is for CDM^HDE. Fig. 2a and 2b are for 
the decay rate of (l)-type. Fig. 2c and 2d for the decay rate of (2)-type, and Fig. 2e and 
2f for the decay rate of (3)-type. 

is shown Fig. 4e and Table 1, there is a good value of Zeq = 1-1 x 10^, which is the 
same order as the observational value z°}' if the interaction between CDM and 5DCDM 
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is included. However, we do not resolve the coincidence problem because there is no 
interaction between HDE and CDM. 

We stress that if one uses 5D CDM x in the brane cosmology instead of radiation, its 
late time evolution is not sizably different from the FRW universe with radiation-matter- 
dark energy. 

Concerning the type of decay rate F, we find that (3)-type is suitable for all interactions 
and thus it could be regarded as the general form. (2)-type works for three cases of 
HDE^CDM, CDM^HDE, and SDCDM^HDE. Finally, (l)-type works for HDE^CDM 
only and it belongs to a very restricted decay rate. 
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Figure 3: (color online) Six graphs for the interaction between HDE and 5DCDM. For 
6^ = 0.2 and c = 1, k = 1 evolution of (green), (red), and (blue) and the 
effective equations of state, u'^ (cyan) and (yellow) with cjj^ = 0. The left column is 
for HDE-^5DCDM and the right one is for 5DCDM-^HDE. Fig. 3a and 3b are for the 
decay rate of (l)-type. Fig. 3c and 3d for the decay rate of (2)-type, and Fig. 3e and 3f 
for the decay rate of (3)-type. 
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Figure 4: (color online) Six graphs for the interaction between CDM and 5DCDM. For 
b"^ = 0.2 and c = 1, k = 1 evolution of (green), (red), and (blue) and the 
effective equations of state, 00]^ (cyan) and (yellow) with uj^ (pink). The left column 
is for CDM^SDCDM and the right one is for SDCDM^CDM. Fig. 4a and 4b are for 
the decay rate of (l)-type. Fig. 4c and 4d for the decay rate of (2)-type, and Fig. 4e and 
4f for the decay rate of (3)-type. 
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